Conductance fluctuations in chaotic systems with tunnel barriers 
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Quantum effects are expected to disappear in the short- wavelength, semiclassical limit. As a mat- 
ter of fact, recent investigations of transport through quantum chaotic systems have demonstrated 
the exponential suppression of the weak localization corrections to the conductance and of the Fano 
factor for shot-noise when the Ehrenfest time te exceeds the electronic dwell time td- On the 
other hand, conductance fluctuations, an effect of quantum coherence, retain their universal value 
in the limit teItd -* 00 , when the system is ideally coupled to external leads. Motivated by this 
intriguing result we investigate conductance fluctuations through quantum chaotic cavities coupled 
to external leads via (tunnel) barriers of arbitrary transparency F. Using the trajectory-based semi- 
classical theory of transport, we find a linear r^-dependence of the conductance variance showing 
a nonmonotonous, sinusoidal behavior as a function of F. Most notably, we find an increase of 
the conductance fluctuations with te, above their universal value, for T < 0.5. Fhese results, con- 
firmed by numerical simulations, show that, contrarily to the common wisdom, effects of quantum 
coherence may increase in the semiclassical limit, under special circumstances. 

PACS numbers: 03.65.Sq, 05.45.Mt, 73.23.Ad 



I. INTRODUCTION 

Since the foundation of quantum physics there has 
been huge interest in the nontrivial transition from the 
quantum to the classical regime. An important obser- 
vation in this context is the Ehrenfest theorem stating 
that the dynamics of quantum mechanical expectation 
values is determined by the classical equations of mo- 
tion [1]. Going beyond expectation values, an Ehrenfest 
time scale has been identified as the time below which 
the quantum time evolution is well approximated by the 
corresponding classical dynamics [2] . The Ehrenfest time 
is the time it takes for the chaotic classical dynamics to 
stretch an initially narrow wave packet to some relevant 
classical length scale such as the system size L. Since the 
stretching is exponential in classically chaotic systems, 
one has 
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(1) 



with the Lyapunov exponent A of the classical dynamics, 
the initial spread h/pF of the wave packet and the Fermi 
momentum pf- 

In recent years, there has been much interest in deter- 
mining the influence of the Ehrenfest time on stationary 
transport quantities such as the conductance [3-8], its 
variance [6, 9] and its behavior under decoherence [3, 10— 
12], shot noise [13-15] and higher moments of the current 
[16, 17], and on time dependent quantities such as the 
spectral form factor [8, 18-20], the survival probability 
[21, 22] and the fidelity [23]. Most of these papers used 
the trajectory-based semiclassical approach to transport 
[24-26] which currently is the method of choice for inves- 
tigating Ehrenfest-time dependences of quantum observ- 
ables [6, 27] . The leading order quantum correction to the 



FIG. 1: A chaotic cavity coupled to external leads via tunnel 
barriers of transparency F < 1 (gray boxes). A classical tra- 
jectory traversing the system is shown by a solid (blue) line. 
The second dashed (red) line on the left indicates a trajectory 
backreflected at the barrier. This occurs with probability 1— r. 



conductance was found in Refs. [3, 5-7] to decay expo- 
nentially with the Ehrenfest time. Qualitatively speaking 
this can be understood by noting first that this contribu- 
tion originates from loop diagrams [see Fig. 3(b)], and as 
such depends on the return probability, and second that, 
neglecting system-dependent nongeneric processes, this 
return probability contribution essentially vanishes for 
times shorter than the Ehrenfest time. Such an intuitive 
interpretation of the Ehrenfest-time dependence does not 
always work: though inherently of nonclassical nature, 
the leading order contribution to the conductance vari- 
ance of systems ideally coupled to external leads turns 
out to be independent of the Ehrenfest time [6, 9]. 

In this paper, we determine semiclassically the 
Ehrenfest-time dependence of the variance var G(E) of 
the conductance for a chaotic system coupled to external 
leads via nonideal contacts modeled by tunnel barriers of 
transparency T < 1. The situation is depicted in Fig. 1. 
The presence of tunnel barriers has the dramatic effect 
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FIG. 2: Dependence of the transmission fluctuations on cavity 
lead coupling F. Fhe function qtb/tc, defined in Eqs. (77) 
and (78), measures the deviation from the universal RMT 
variance. In this graph the parameter a is plotted as a func- 
tion of the tunnel barrier transparency T for time-reversal 
symmetric systems. The solid red line gives the analytical 
prediction, Eq. (78), and the blue dots are results of numer- 
ical simulations. Error bars indicate the standard deviation 
over the ensemble of calculated data. 

that varG(-E) increases or decreases with te, depend- 
ing on the value of T. For T < 0.5, we even observe an 
enhancement of the variance above the universal value 
in the presence of time reversal symmetry for equal lead 
widths of var(G) RMT = (l + (l-r) 2 )/8 [28] upon increas- 
ing te- This is very surprising, given the quantal nature 
of the conductance fluctuations. In the range 0.5 < T < 1 
we find a reduction of vsltG(E) as te increases that is 
strongest around T ~ 0.8, and recover the T^-independent 
behavior of var G(E) at T = 1 observed in Refs. [6, 9]. 
The precise dependence on T is depicted in Fig. 2. 

To obtain these results we first, in Sec. II, introduce 
the semiclassical approximation to the conductance vari- 
ance. Then, in Sec. Ill, we analytically calculate the 
T-dependence of the diagrammatic contributions to the 
variance, to leading order in the total number of open 
channels and linear in the Ehrenfest time. We list all rel- 
evant diagrams and calculate their contributions. Given 
their number, we also identify the most relevant ones and 
specify the range in T where they are particularly impor- 
tant. In Sec. IV we present numerical results that confirm 
our analytical results and finally conclude in Sec. V. 



II. SEMICLASSICAL APPROACH 

Within the scattering approach to quantum transport 
[29, 30], the energy averaged conductance G(E) (in units 
of 2e 2 /h) can be expressed in terms of the transmission 
t as 

G(£) = (Tr(ttt)} (2) 



with (••■) denoting an average over an energy window that 
is classically small but quantum mechanically large. This 
leads to the following expression for the variance 

vaxG(E) = ([Tr(tt t )] 2 )-(Tr(tt t )) 2 . (3) 

The scattering matrix elements are related by the Fisher- 
Lee relation [31] to the projection of the Green function 
onto the transverse directions in the leads. Performing 
the projection to leading order in h and approximating 
the Green function semiclassically, one obtains 

ta,*--^ E (4) 

V 1 H 7 ( a ^fc) 

with the Heisenberg time Th, the time conjugate to the 
mean level spacing. Here the sum is over the scattering 
trajectories 7 which connect channel a in the entrance (or 
say left) lead and channel b in the exit (or right) lead in 
Fig. 1. The summands contain rapidly oscillating phases 
depending on the classical actions 5 7 of the considered 
classical trajectories, and classical stability prefactors A 7 
whose precise form is given for example in [26]. 

Inserting Eq. (4) into (2) we obtain the semiclassical 
expression for the conductance, 

G(E)»l±Z E A 7 A;,e^~ s A. (5) 

\ 1 H a,fc 7 , 7 '(a^h) / 

Using (4) in (3) yields the semiclassical expression its 
variance, 

var G(E) « /-L^ E A^A^, (6) 

\ 1 H a,b 7j7 '( a ->b) 

xe i(S7-Sy+.s«-%)J 

~(^E E A 7 A;,e^- s A 2 

\ 1 H a ,b 7,7'(o->6) / 

with the channel sums in Eqs. (5) and (6) running over all 
open lead channels (JVl in the left and JVr in the right 
lead). If we consider contributions in the first term in 
Eq. (6) where 7 and 7' form a correlated pair (with self- 
encounters) and £ and £' form a separate correlated pair, 
we simply recreate the second term. We can thus remove 
the second term in the above equation by removing such 
pairs from the semiclassical treatment of the first term. 
In terms of trajectories we then obtain 

varG(£) * ((-^E E A,A;,A^Al 

xe7 Hs,-s T , + s,-s e ,)\\ (7) 



where the trajectories 7, 7' go from channel a in the en- 
trance lead to channel b in the exit lead. Likewise trajec- 
tories go from channel c to channel d. Because we 
have removed terms from correlated trajectories where 
7 fa 7' and £ rs £' [this restriction is denoted by the dou- 
ble bracket in Eq. (7)] we are left with quadruplets where 
all four trajectories interact through encounters. 

Before performing the energy average the approxima- 
tions for G(E) and var G(E) in Eqs. (5) and (7) are 
rapidly fluctuating as a function of energy for h -*■ 0. 
Thus only contributions from very similar trajectories 
survive the average. The classical contribution to Eq. 
(5) results from equal trajectories 7 = 7', the so called 
diagonal approximation [24, 26], for an illustration see 
Fig. 3(a). Here Eq. (5) yields 



G(E) 



[3(a)] 



1 



1 H a ,b 7 ( a ->o) 



(8) 



From here on, superscripts refer to the corresponding fig- 
ure. To evaluate the remaining 7-summation the sum 
rule for open systems [25] that transforms the sum over 
orbits into an integral over their durations is applied 
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dte~ t/TD . 



(9) 



Here td is the dwell time of the cavity, i.e. the typical 
time a classical particle remains inside the chaotic sys- 
tem. This is given by t d = T H /N with N = N h + N R . 
This finally yields 



G(E) 



[3(a)] M 



Nt, + Nt 



R 



(10) 



Quantum corrections to this result are obtained from 
pairs of slightly different trajectories. The considered 
trajectories are almost identical differing only in how 
they are connected within self encounters, for the pair 
considered in [25], see Fig. 3(b). Here the orbits pos- 
sess close self encounters with one orbit crossing and the 
other anticrossing there leading to a different direction 
of traversal of the closed loop. Considering this pair in 
Eq. (5) leads to the leading order quantum correction to 
the conductance. To determine its contribution the ac- 
tion difference between the partner trajectories and the 
number of crossings needs to be determined. The calcu- 
lation is done here within the phase-space approach, in 
the context of the conductance it was first performed in 
Ref. [32]. We will follow the latter approach throughout 
this paper. There a Poincare surface of section is consid- 
ered inside the encounter region and the difference along 
the stable and unstable directions of the piercing points 
of the two stretches, s and u, respectively, is used to 
characterize an encounter. In terms of these coordinates 
the action difference for the orbit pair in Fig. 3 is given 
by AS = su [32]. The weight w(s,u) that additionally 
depends on the duration of the orbit T measuring the 
number of encounters is obtained from the ergodicity of 




FIG. 3: Schematic drawing of trajectory pairs contributing 
semiclassically to the energy averaged conductance (a) A pair 
of identical trajectories leading to the classical contribution 
to the conductance, (b) The pair depicted here differs at 
a self-encounter and leads to the weak-localisation quantum 
correction to the conductance, (c) An additional pair of orbits 
that needs to be considered in the presence of tunnel barriers. 
This diagram is obtained from (b) by shrinking the left link 
until the encounter touches the tunnel barrier. An analogous 
configuration can be formed by removing the right link from 
Fig. 3(b). 



the flow as [32] 



w(s, u) 



(T 2t enc ) 

2o2ionr 



(11) 



where t cnc = 1/ A In (c 2 /| sit] ) is the duration of the en- 
counter. In general, the action difference AS and the 
weight w(s,u) depend on the trajectory configuration 
considered. In total we obtain for the quantum correc- 
tion SG(E) resulting from the diagram in Fig. 3(b) 



8G(E) 



[3(b)] 



f° f°° iAS ( T ~«e 

/ dsdu dTw(s,u)e h e T o 
J-c Jit, 



(N L + N R y 



(12) 
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with te = l/Aln(c 2 /ft). In the first line additionally the 
survival probability correction during the encounter [32] 
is taken into account. During the encounter the stretches 
are so close that the orbit escapes either during the first 
stretch or does not escape at all leading to the enhanced 
survival probability in Eq. (12). The s, it-integrals in Eq. 
(12) are performed as described in [6, 8], by substitut- 
ing su = c 2 x and a = c/u. The cr-integral then essen- 
tially cancels the i cnCi in the denominator and the x- 
integral yields (after a partial integration) the contribu- 
tion -l/(T D T„)e- TE/TD . 

To treat orbits differing at several places, Ref. [32] in- 
troduces the splitting of the orbit into encounters and 
links. During the encounters the orbits are close to them- 
selves but the orbit and its partner are differently con- 
nected. Due to the exponential separation of neighboring 
trajectories in the chaotic case these last essentially an 
Ehrenfest time, as will become clear from the calculation 
below. The stretches are connected by the links, where 
'links' denote the long parts of the trajectory where the 
trajectory and its partner are essentially identical (up to 
time reversal symmetry) . With a suitable change of vari- 
ables in the calculation, in the RMT limit te -*■ 0, one 
can treat different encounters as distinct and separate the 
semiclassical contribution into a product of contributions 
over each of the links and encounters. The total contri- 
bution can therefore be obtained by diagrammatic rules 
[16, 32]. 

Away from this limit, and for the Ehrenfest-time de- 
pendence, the encounters may start to overlap and for 
the conductance variance the trajectories can be seen to 
meet and surround periodic orbits trapped inside the sys- 
tem, i.e. they have encounters with periodic orbits, and 
must be treated as part of a continuous family [6]. 

When we include tunnel barriers (as in Fig. 1), three 
main changes occur that were originally described in [34] : 

1. The particles enter and leave the cavity with the 
probability T, which leads to a factor T 2 for each 
trajectory pair. 

2. While without tunnel barriers every particle that 
hits the lead leaves the system, now only the ratio 
r of the particles hitting the lead leaves the cav- 
ity. For the links, the effective dwell time is there- 
fore ro/r and the dwell time in the exponential in 
(9) should be replaced by this effective dwell time. 
However, if trajectory stretches are correlated, as 
they are during encounters, then the whole config- 
uration is lost if just one stretch of the encounter 
leaves the system. This happens with a probability 
p n = 1 - (1 - r) n for n correlated stretches. The 
dwell time in such a situation is therefore replaced 
byT D /[l-(l-T) n ]. 

3. Additional encounter diagrams become possible; 
for an example which contributes to the energy 
averaged conductance G(E), see Fig. 3(c). In 
this case one encounter stretch can be moved into 




FIG. 4: (a) A diagram possessing two 2-encounters in a row. 
This diagram does not require time-reversal symmetry unlike 
the corresponding diagram in (b). 

the lead forcing the other to be backreflected at 
the opening. Note that configurations where both 
stretches are backreflected at the opening are al- 
ready taken into account by the modified dwell time 
explained above. 

Although the effective dwell times are altered by the 
tunnel barriers, the action difference and weight func- 
tions are unaffected so that in the RMT limit (te -»■ 
0) contributions can still be obtained by diagrammatic 
rules: The contribution of each link is now given by 
(TN)^ 1 . The stretches of an encounter of n orbits yield 
—p n N. For the Ehrenfest-time dependence, however, 
these changes render the calculation of the contributions 
to the conductance variance much more difficult com- 
pared to r = 1 : Due to the discontinuous form of the ef- 
fective dwell time, the contributions from diagrams with 
a different number of surroundings of trapped periodic 
orbits need to be split and treated separately. Also the 
possibility for encounter stretches to be backreflected at 
the tunnel barriers increases the number of diagrams con- 
siderably. 

III. DIAGRAMMATIC CONTRIBUTIONS 

Here we calculate the leading order in N contributions 
to the variance of the conductance for non-zero Ehrenfest 
time in the presence of tunnel barriers. We show all the 
relevant diagrams and calculate their contributions. The 
results given here are valid in the unitary case; results 
for the orthogonal case can be obtained by multiplying 
the total by a factor of 2. 

A. Discrete encounters 

We start with the contributions important in the 
RMT limit which allow us to recover the RMT re- 
sult. The corresponding RMT calculation was performed 
in [28] by Brouwer and Beenakker. First we consider 
two 2-encounters (encounters involving two trajectory 
stretches) in a row, see Fig. 4(a). This diagram also oc- 
curred in the T = 1 treatment of [6], but we will explain 
how for r + 1 other diagrams with backreflected stretches 
can be derived from this one. 

When the encounters are inside the cavity and not 
touching the tunnel barriers, as depicted in Fig. 4(a), we 
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start with Eq. (7) and use the sum rule (9), the action 
difference AS(s,u) = S\U\ + s 2 u 2 (where the subscripts 
refer to the different encounters) and the weight function 

[32] 



j(s,u) 



(Tl - tcnc,l _ Wc,2) (T 2 - tcnc,l ~ t cnCj2 ) 
4f2 2 £ cnc i t cnc 2 



(13) 

containing the durations of the two trajectories indicated 
by solid lines in Fig. 4(a), T x and T 2 , respectively and the 
durations of the two encounters of the trajectories [33] 



te 



i€{l,2} 



(14) 



with a classical constant c again of the order one. From 
Eqs. (6) and (9) we write 



var 



G [4] 



4 /\f2 ]\t2 6 



H i=l J0 



dt, e- ru/TD / dsdu 

J-c 



3 -p2(tcnc,l+tcnc,2)/T D 



^enc.l^ 



1 

Q 2 

(15) 



cnc, 1 ^cnc,2 



with the superscript again referring to the corresponding 
figure: Fig. 4. As explained in Sec. II, the trajectory 
quadruplet leads to the overall factor T 4 while the sum 
over possible channels provides N^N^. The six links have 
an effective dwell time of td/F, while each 2-encounter 
experiences the dwell time td jvi as explained at the end 
of the last section. The s, u-integrals are performed, as 
described after Eq. (12), by substituting SiUi = c 2 Xi and 
(7j = c/ui. Each Xi-integral yields (after a partial integra- 
tion) the contribution -p 2 / (tdTh) e~ P2TE l TD as already 
obtained for the conductance for T = 1. Finally the U- 
integrals yield 



var 



G [4] = 



NlN* R 



(N L + N R ) 4 ^ 



2 

P2 r -2p 2 T E /T D 



(16) 



which generalizes the result for T = 1 from [6] . The repre- 
sentation of Eq. (15), and the integrals to arrive at (16), 
also nicely illustrate how the diagrammatic rules intro- 
duced above arise in this context. 

We now turn to the new diagrams that arise due to 
the tunnel barriers where some of the links are shrunk 
until an encounter touches a barrier. First we consider 
the case in which just one link connecting the encounter 
to the opening is removed and the corresponding en- 
counter stretch now starts in the opening; as in the 
example in Fig. 3(b). As the stretches during an en- 
counter lie very close to each other, the other encounter 
stretch has to be backreflected at the opening (so only 
one link is lost). This contribution can therefore only 
exist for V + 1. The changes in the analytical calculation 
mainly affect the weight function, see also [8, 21, 22, 34]. 
If we shrink a link on trajectory 1 then, compared to 
Eq. (13), the orbit of duration T\ now only involves two 
links so the factor (7\ - t en c,i _ Wc,2) /2 is replaced by 



(Ti - £ e nc,i ~ Wc.2)- Moreover, for the encounter which 
touches the barrier, t cnc ^ is replaced by an integration 
variable t' that is integrated from zero to the t enc _i de- 
fined in (14). This variable measures the length of the 
encounter that remains inside the system, i.e. which has 
not yet been moved into the lead. Performing again the 
steps described after (14) yields an expression similar to 
(15) but with 5 instead of 6 link factors and t C nc,i in the 
exponential in the second line replaced by an integration 
variable t' that is integrated from zero to £ cnC; j. Because 
half of the encounter is backreflected at the tunnel bar- 
rier, additionally this contribution is multiplied by (1-r). 
The contribution from the lower limit of the t'-integral is 
zero [8, 17], leading to 



var 



G [4-ii] = 4NlNj (1-I>2 c 

(^Vl + ^Vr) 4 r 



(17) 



with the '-11' denoting that one link was removed. The 
prefactor 4 is due to the fact that there are four such links 
we can remove. Analogously, the contribution where two 
links that connect two different encounters to the opening 
are removed is 



var 



G [4-21] 



4A T 2 iV, 



^^ T (i_r)V 2p2rE/rD , (is) 



(A^l + JV R ) 



where there are again four possibilities of picking two 
links to remove. 

Additionally we can, and this is a possibility which 
also exists for T = 1, remove both links connecting the 
same encounter to the opening. This means that the 
trajectories tunnel straight into the encounter. As the 
encountering orbits are so close together this means that 
a = c or b = d so that there is only one channel summa- 
tion in the lead where the encounter touches. Also we 
have two links fewer and one integral over the part of the 
encounter that remains inside the system, yielding 



varG [4-21( S )] = _ 



(NlN R + N L N p 
(N L + N R f 



-p 2 e 



-2p 2 T E /T D 



(19) 



where the additional ' (s) ' in the superscript indicates that 
two links were removed at the same encounter. We can 
further remove one link from the other encounter to ob- 
tain 



var 



r , , n 2 (N?N r + JVlJVh ) „ . 

G [4-31] _ V L, ti ^ rci y)c P2Te ^ Td 

(N L + N R f 



(20) 



Finally, when all four links connecting the encounter to 
the leads are removed, we have 



var 



G [4-4i] = N h N R t 2 c -2 P2 t e /t d 
(A L + A R ) 2 



(21) 



With time reversal symmetry, however, we can also 
have Fig. 4(b) where, because channels a and c are in the 
left lead and channels b and d in the right lead, we cannot 
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FIG. 5: (a) A diagram possessing two independent 2- 
encounters with a periodic orbit (dashed-dotted line). While 
this diagram does not require time reversal symmetry, the 
corresponding diagram (b) does. 




FIG. 6: (a) A diagram with one 3-encounter with a periodic 
orbit (dashed-dotted line). Diagram (a) does not require time- 
reversal symmetry while diagram (b) does. 

shrink more than two links simultaneously. Similarly, we 
cannot remove two links from the same encounter. Since 
P2 = 2r - r 2 , the contributions (19-21) actually cancel 
so that the diagram in Fig. 4(b) and the ones obtained 
by cutting links provide the same contribution as the di- 
agram in Fig. 4(a). Time reversal symmetry then still 
gives a factor 2 in this case, while for all the following 
cases, diagrams related by time reversal symmetry pro- 
vide the same contributions directly. 

Next we treat the diagrams with two 2-encounters 
which lie along a trapped periodic orbit; depicted in Fig. 
5. Although the links connect the encounter stretches in 
a different way compared to Fig. 4 these diagrams again 
contain two 2-encounters and the same number of links. 
The possibilities for shrinking links are identical as for 
Fig. 4(b) and so each diagram provides the same contri- 
butions as above and the same total contribution as the 
configuration in Fig. 4(a). They are multiplied, however, 
by a factor 2 since the orbits here have two possibilities 
to traverse the enclosed periodic orbits (schematically we 
can also reflect the diagrams horizontally). 

The last relevant diagram type to obtain the RMT 
result is shown in Fig. 6. This diagram contains one 
3-encounter with a periodic orbit. Its contribution is cal- 
culated in an analogous way to Eq. (15). Here we have 
one link fewer, one encounter instead of two and involv- 
ing three rather than two orbit stretches. Therefore we 
have 

varG [6] = - 2N * N * gyw^, (22 ) 

(N L + N R ) 4 T 



Alternatively this result can be obtained from the con- 
tribution of a 3-encounter, K\, in the Appendix of [8] by 
including the correct dwell times as well as the contribu- 
tions from the links. 

Also in this case encounter stretches can be shrunk 
and removed: First removing one link connecting the 
encounter to the opening yields 

varG [6-iq = ^M_ (1 _ r) 2 e - P3 r E /r D (23) 

and second removing both links: 

rfi „ n 2 (N h N2 + N?N K ) , 
varG [6 - 2i] = V R V^r(l-r)e- p3TE/TD . (24) 

(N L + N R y 

Having calculated all the contributions in the RMT 
limit, we can obtain the RMT result by setting te = 0. 
When summing the the results in Eqs. (16-24), we obtain 

rmt = n l n r t 6 j jvgjvg (4 -sr + er 2 ) 

+ (N[N R + N h N%) (2 - 2r + r 2 ) 

+ K + ^ L 4 )(2r-2r 2 )], (25) 

which agrees with the RMT prediction in [28]. In Ap- 
pendix B, we use these diagrams to obtain the RMT 
result when each channel has a different tunneling prob- 
ability. 

B. Periodic orbit encounters 

Having gone through all the diagrams that contribute 
at zero Ehrenfest time, we now turn to those diagrams 
whose contribution vanishes at zero Ehrenfest time. For 
these contributions, the periodic orbits in Figs. 5 and 6 
become important and we now view those diagrams as 
trajectories which have an encounter with the periodic 
orbit, rather than with each other. For example, in Fig. 
6 we could start with the solid trajectory which passes 
from channel c to d and the dashed trajectory from a to 
b and build the rest of the diagram from those starting 
points and the periodic orbit. Both of those trajectories 
encounter the periodic orbit once. In the semiclassical 
treatment of Fig. 6 above, it was implicitly assumed that 
these encounters occur at the same point along the pe- 
riodic orbit. The resulting 3-encounter can therefore be 
considered as an 'aligned' 3-encounter, but for the fur- 
ther Ehrenfest time dependence we also need to consider 
the situation where the two encounters with the periodic 
orbit occur at different points along the periodic orbit 
but still overlap. In this case we have a 'non-aligned' 3- 
encounter, while when the encounters no longer overlap 
we return to the two separate 2-encounters of Fig. 5. 

We now derive the Ehrenfest-time dependence of a 
non-aligned 3-encounter, whose base trajectories are de- 
picted in Fig. 7(a). In order to obtain the complete tra- 
jectory quadruplet that contributes in (7) we first need 
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(a)-... 




FIG. 7: Periodic orbit encounters which contribute only for 
nonzero Ehrenfest time, (a) the base trajectories for a non- 
aligned 3-encounter and (b) for encounters which overlap at 
both ends thus enclosing the periodic orbit (depicted dashed 
dotted). The encounter stretches are shown thick (gray), 
while the links connecting the encounter stretches to the open- 
ing are indicated by solid (blue) and dashed (red) lines. To 
obtain the complete quadruplet of trajectories (with a small 
action difference), an additional traversal of the periodic or- 
bit must be included with one of the base trajectories and 
included with the other base trajectory for the partner tra- 
jectories. 



to include an additional traversal of the periodic orbit 
with one of the base trajectories to obtain the original 
trajectory pair [which has positive action in (7)]. The 
partner pair (with negative action) is then created by in- 
cluding the extra traversal of the periodic orbit with the 
other base trajectory. In this way we recover a diagram 
like Fig. 6(a) from Fig. 7(a) and a small action difference 
(the action of the periodic orbit itself cancels). After 
writing the contribution of a non-aligned 3-encounter in 
an analogous way as for the aligned 3-encounter, as ex- 
plained in Eqs. (13-15), as product of link and encounter 
contributions, it can be evaluated by making use of the 
results for K 2 in the Appendix of [8]. Of course, with 
the tunnel barriers, the dwell times must be modified 
compared to [8]: The link dwell time is modified by the 
factor 1/r, while the parts of the encounter where only 
one trajectory is correlated to the periodic orbit — these 
contributions are called fringes in [8] — have the fac- 
tor l/p2 since we have two stretches close to each other. 
Likewise, when both trajectories are correlated with the 
periodic orbit we have three stretches in total and the 
corresponding factor With these corrections, the 

contribution is 



• G [7(o)] 



(N L + N R yr(2 P2 -p 3 ) 



(< 



-P3Te/Td 



-2p2 t e /t d 



(26) 



In order to proceed to the additional diagrams that 
arise from touching the tunnel barriers, we first recon- 
sider this contribution in detail along the lines of [17]. 
As for Eqs. (13-15), we start with the weight function 
for two base trajectories encountering a periodic orbit, 



as explained in [6, 8] 



d*l / d *2 7^ 

Jo 



T 2 -t e 



to 



-. f dr p f dt'. 



(27) 



Here t\ and t 2 are the durations of the links of the base 
trajectories that connect the periodic orbit encounters to 
the lead, while t p is the period of the periodic orbit. The 
integral over the period corresponds to summing over all 
periodic orbits which can be encountered. The encounter 
times are as in (14), but using the stable and unstable 
distances between the encounter stretches and the peri- 
odic orbit itself. Finally, t' measures the time difference 
between the midpoints of the two encounter stretches in 
Fig. 7, and the i'-integral covers the different arrange- 
ments of the stretches relative to each other. The limits 
of the two last integrals in (27) are not given as they de- 
pend on the specific configuration considered below. The 
first two integrals in (27) can again be transformed, to- 
gether with the integrals from the sum rule (9), into a 
product of link and encounter contributions. The gen- 
eral expression for the contributions to the variance from 
diagrams containing an enclosed periodic orbit that does 
not touch the lead becomes 



varG po 



T*N?N? ! 



t 



i^nx 

J dr p J di'e-^ 7- "*-"' 1 '*-"'* 



enc,l^enc,2 



(28) 



Like the limits of the t p - and i'-integrals, the function P 
determining the survival probability along the periodic 
orbit (and the encountering trajectory stretches) will be 
specified for each contribution separately below. 

For example, for the non-aligned 3-encounter, we can 
arrange the different alignments in terms of the durations 



of the two encounters. We let t P 



denote the longer 



encounter and i cnc ,min the shorter. In Fig. 7(a), we do 
not yet allow the encounter stretches to surround the pe- 
riodic orbit so we impose the restriction t p > f CM , max . 
Then we can first consider the case that the shorter en- 
counter lies inside the longer. We will also refer to this 
case later as a 'generalized' 3-encounter. The time differ- 
ence between the midpoints of the encounters therefore 
satisfies \t'\ < (£ C nc,max _ ^cnc,min) /2, while the survival 
probability function in (28) is given by 



P 



T(t p - 1, 



p ^enc.max 



) , P2 (ft 



enc.max ^enc^mirw 



TD 
P3^cnc,min 
TD 



TD 



(29) 



The different terms simply correspond to the parts of the 
periodic orbit which are followed by one, two and three 
trajectory stretches respectively. As this survival prob- 
ability is independent of t', the ^'-integral simply yields 
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(ienc,max - We, min)- Performing the remaining integrals, 
we obtain 



2N£N£ P3 -2 P2 c _ P3Te/td 

(n l + iv R ) 4 r 



(30) 



Of course, the range of t' includes 0, the case in which 
the encounters are perfectly aligned. Eq. (30) therefore 
includes the contribution of the 3-encounter in Fig. 6(a), 
and for just the non-aligned 3-encounter we subtract (22) 
to give 



(31) 



Now we consider the case in which the shorter en- 
counter is no longer fully inside the longer, but where 
the two encounters still overlap. This can be further 
separated according to the total length of the encoun- 
ters, t oncs = i C nc,i + Wc,2, compared to the periodic orbit. 
When t p > tones, we have the range 



1^ I ^ [(^cnc,n 



-t, 



cnc,min 



) /2, ^encs/2] • 



(32) 



and the survival probability 



Tt p (P2-T)(t 



p = __L + ^ -' I + \£ 

T D T D \ 2 



(P3-P2) (t 



(^p+l^'l) 



(33) 



where the terms are now expressed as corrections due to 
additional correlated stretches. Integrating the exponen- 
tial depending on this survival probability in (28) with 
respect to t' yields 



2t d 



p 3 - 2p 2 + r 



-[(p2-r)i encB +rTp]/r 



g — [p3^cnc,min +P2 (^cnc,max ^enc , min ) + T (Tp icne , max )]/t£> 

Performing the remaining integrals in (28) leads to 
4AT£iV£ 1 



(34) 
)■ 



(jv l + n r ) 4 r (p 3 - 2p 2 + T) 
x [p 2 2 e- 2p2TE/TD -p 2 (P3-P2 - r) c - (p3+r)TE/TB ] 



(35) 



If the encounters are longer than the periodic orbit, 

^enc.max 

< t p < tones, but we still do not allow them to 
overlap at both ends, we have the restriction 

|t'| 6 [(t 

cnc,max ^cnc,min) /2, Tp ~ toncs/2] , (36) 

while the survival probability remains as in (33). The 
t'-integral yields 



2t t 



P3 -2p 2 + r 



- [ (P3 ~P2 ) tencs + (2p 2 ~P3 )t p ]/td 



(37) 



g [P3^cnc,min +P2 (^cnc,max ^enc , min ) + T (Tp ^enc , max ) ]/t£> ^ 



which finally leads to 



4N l N R 



1 



(JV L + JV R ) 4 i(2p 2 -p 3 )fe 



X[(P3-P2)P2C- P3TE,TD -pjc 



2 P2 + r) 

2p2T E /r D J 
[(P3-P2)P2C- P3TE/TD 



T(p 3 -2p 2 + T) 

~P2 (p3-p 2 +r)c- (P3+r)TE/T °]}. 



(38) 



As can be easily checked, the sum of the contributions 
of the non-aligned 3-encounter in Eqs. (31,35,38) equals 
the contribution calculated directly in (26). The reason 
why we have discussed this more complicated route, is 
that we can use it to easily calculate the contributions 
when we start to shrink links and allow the encounter to 
touch the tunnel barriers. For example, if we remove one 
link, then by performing analogous steps as explained 
before (17), we obtain 



varG [7(a)-ll] 



% N t N R 



p 2 (i-r) 



[i + (i-r)] 



{Nl + Nl) (2p 2 -p 3 ) 



where in the square brackets, the 1 results from config- 
urations where the enclosed periodic orbit touches the 
tunnel barrier where only one encounter stretch is corre- 
lated with the orbit, and the (1 - T) from configurations 
where the periodic orbit touches the tunnel barrier while 
both encounter stretches are correlated with it. Likewise, 
if two links connecting different stretches to the lead are 
removed, we obtain the contribution 



varG [7(a)-21] = 



4NlN R 



r(i-r) 2 (2-r) 2 



x 1 e 



- P 3T E /t d 



(2p 2 -p3) 
-2 P2 tb/td\ (40) 



Next we turn to the configuration in Fig. 7(b) with 
encounter stretches overlapping at both ends. We still 
have t cnc max < T p < tones, but a different restriction on t'\ 



|t I 6 \jp tcncs/2, 7"p/2, ] ■ 



(41) 



The survival probability is again independent of t', and 
given by 



P = 



(P3 ~P2) t, 
TD 



■ ncs (2ff2 



■P3)T P 



TD 



(42) 



so that the t'-integral yields (t e 
remaining integrals we obtain 



- T p ). Performing the 



var 



G [7(b)] 



2N l N R 

(N L + N R ) 4 



Pi 



(2p 2 -p 3 ) 

-2p 2 T E /T D _ e -p3T E /T D ^ 

P2 (P3-P2)t E „- P3Te / Ti , 



x 1 e 



(2p 2 -Ps)t d 



(43) 
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Again the encounter stretches can be brought to the lead 
by first removing one link 



varG [7(b)-ll] 



p 2 r(i-r) ; 



(A L + A R ) 4 L (2p 2 -p 3 ) 2 

P3 r(i-r) 2 

2(2p 2 -p 3 )T D 



(44) 



and second by removing two links connecting different 
stretches to the opening 



varG [7(b)-21] 



8N l N R 

(A L + N R ) 4 

, r 2 (i-r)S c -^ /Tn 

(2p 2 ~P3)t d 



r 2 (i-r) 4 
(2p 2 -p 3 y 



) 



(45) 



C. Touching both leads 

As the durations of the encounters are of the order of 
the Ehrenfest time, for vanishing Ehrcnfcst time we only 
considered the situation where the encounter stretches 
could be partially reflected from the tunnel barriers in 
one lead. However, for increasing Ehrenfest time config- 
urations where encounter stretches are partially reflected 
from the tunnel barriers in both leads, i.e. where they 
touch the opening at both ends, become important. If 
one of the 2-encounters in Fig. 4 were to be partially re- 
flected at both ends then one of the links between the two 
encounters would need to tunnel through the barrier and 
exit the system so that the rest of the diagram could not 
be completed. In Figs. 5 and 6, however, as long as the 
trajectory stretches which follow the periodic orbit are 
reflected at the tunnel barriers and remain in the sys- 
tem, we can allow the other links to tunnel through the 
barrier and exit the system. For the base trajectories in 
Fig. 7(a) this means we can allow both links of one of the 
base trajectories, on either side of the same encounter, 
to shrink into the start and end lead. We start with the 
configuration in Fig. 5 and for the calculation we note 
that the contributions of the different orbit parts in (7) 
are multiplicative [32]. We can therefore reconnect the 
orbits in such a way that they split into parts whose con- 
tributions have previously been calculated. First there is 
the remaining base trajectory, with a 2-encounter with 
the enclosed periodic orbit and two links connecting it to 
the opening. This contributes the factor 



A L A R 



( A L + A R ) 



TP2G 



(46) 



The rest of the diagram involves the periodic orbit it- 
self and the encountering stretch that tunnels through 
to start and end in the leads. The contribution is calcu- 
lated by summing over all enclosed periodic orbits, using 



the sum rule (9). Remember, that when allowing the en- 
counter to move into the lead an additional time integral 
occurred measuring the part of the stretch that lies still 
inside the system. Now there are two time integrations 
representing the amount of the encounter which is cut 
short in each of the two leads. The first time integral 
cancels t cnC:i , while the second essentially yields a factor 
[1 - exp (-P2t enc ,i/TD)]- For the details of the calcula- 
tion we refer to [17, 27]. This part of the diagram then 
contributes 



A L A R 



(n l + N R y P2 



(!_£)! (l-e-w/ri,) 



(47) 



so that this configuration of Fig. 5 altogether yields 

varG [5-21 (s) ] = _ 

x(l- c - p2TE/TD ), (48) 



4N t N R , x _ r x2 e -p 2 r B /r D 

(A L + A R ) 4V 



where the '(s)' in the superscript denotes that the two 
links were removed along the same base trajectory from 
the same encounter stretch. One factor 2 in the last 
equation derives again from the mirror symmetry of this 
configuration explained after Eq. (21), the other from the 
fact that each of the two stretches can touch the opening 
at both ends. 

Furthermore, one of the links of the other base tra- 
jectory may also be shrunk so that the other encounter 
tunnels into the lead at one end. This contribution is 



var 



G [5-31] 



8A 2 A 2 



(A L + Ar) 4 P2 

x(l-e~ P2TE/TD ). 



r(i-r) 3 P2Te/td 



(49) 



Removing all 4 links, so that both encounter stretches on 
the base trajectories tunnel into the leads at both ends, 
is also possible. The contribution is simply given by the 
square of (47). 

We can repeat this treatment for the aligned 3- 
encounter of Fig. 6. If one encounter stretch tunnels into 
the leads at both ends we obtain 



varG [6-21 (s) ] 



8N t N R 



(A L + A R ) 4 P3 



F(1 r)4 (l- c -^/^), 



(50) 



where, because of the alignment and proximity of the two 
encounter stretches, the periodic orbit and the other en- 
counter stretch must be backreflected at the tunnel bar- 
riers. Allowing the ends of the second encounter stretch 
to progressively tunnel through into the leads as well, we 
have, if it tunnels into the lead at one end 



var 



G [6-3i] 



4(a l a 2 +a 2 a r ) r 2 (i-r) 3 

(A L + A R ) 3 P3 
x(l-c- p3TE/TD ) 



(51) 
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be expressed as 



FIG. 8: A periodic orbit encounter with fringe correlations, 
i.e. correlations between the two trajectories that encounter 
the periodic orbit (shown dashed dotted) with each other, 
but not with the periodic orbit itself. The fringe regions are 
indicated by the black lines perpendicular to the trajectories. 
The encounter stretches are shown thick (gray), while the 
links connecting the encounter stretches to the opening are 
indicated by solid (blue) and dashed (red) lines. 



and at both ends 

2N L N n r 3 (i-r) 2 



, £t[6-41] 



(N h + N R y 



P3 



(52) 



Such configurations can be also considered for the dia- 
grams in Fig. 7 and their contributions can be calculated 
analogously. To simplify the results, however, we will 
later perform an expansion of the contributions in pow- 
ers of the Ehrenfest time and only retain terms up to 
linear order. As the results for the diagrams of Fig. 7 are 
all of higher order in te than the linear one, we will not 
focus on their explicit form here. 



D. Encounter fringes 

A further effect, and one that actually causes the in- 
dependence of the conductance variance of the Ehrenfest 
time for T = 1, are correlations during fringes near pe- 
riodic orbits. Encounter fringes refer to regions where 
the two base trajectories which encounter the periodic 
orbit are correlated with each other, but are no longer 
correlated with the periodic orbit itself; see Fig. 8 for a 
schematic depiction. 

As the two base trajectories leave the periodic orbit 
correlated with each other, we can consider the encoun- 
ters with the periodic orbit to be aligned. The encounter 
of both trajectories with the periodic orbit has length 
t Gnc while the fringes, during which the base trajectories 
are correlated have lengths t s and t u where the subscripts 
refer to the fact that the stable and unstable distances 
between the base trajectories must be small for them to 
remain correlated. Here, we treat the case in which the 
encounter length is shorter than the periodic orbit. By 
this we generalize the calculation of [6] to T + 1. There, 
it was shown that the contribution derived from (28) can 



var 



j / ds'du' / 



(*hy(N h + N K ) 



ds du / din 



[(P3-r)te„c+rr p ]/r D 




x cos 



(53) 



where b is again a classical constant of order unity. The 
stable and unstable coordinates are defined as in [6] : As 
correlation effects away from the periodic orbits get only 
important when the encountering orbits approach and 
leave this at classically close points one first characterize 
the correlation of all three orbit pieces (i.e. the periodic 
orbit and the two encountering orbits) by the coordinate 
u. Furthermore the difference between the stable and the 
unstable coordinates of the encountering orbits at the be- 
ginning and the end of the encounter is denoted by s' and 
u' . These coordinates characterize the time the encoun- 
tering orbits remain correlated before and after they are 
correlated with the periodic orbit, respectively. In (53), 
we included the correct dwell times for the different parts 
of the trajectory: Tt>/r for the isolated links, Tu/p2 for 
the fringes and Tjj/p^ for the encounter stretches corre- 
lated with the periodic orbit. Following the steps in [6], 
this contribution evaluates to 



var 



2N£N* 



Pi 



(N L + N B ,fT(2p.. 
x ( e - p3TE/TD 



2-P3) 

-2p 2 T E /T D ^ 



(54) 



Without tunnel barriers, as shown in [6], further contri- 
butions can be obtained when the fringes start in the 
leads. Similarly, with tunnel barriers, we can have the 
fringes tunnel into the leads. This leads to an additional 
time integral over the duration of the fringe that remains 
inside the system. When one fringe tunnels into the lead, 
two links are removed and the contribution becomes 

varG [8-2i( s) ] = HNlN R+ N L N*) T P2 

x ( e -P3T E /r D _ e -2p 2 T E /r D J _ ( 55 ) 

The '(s)' in the superscript refers to the fact that the 
links were removed from the same side of the encounter. 
When both fringes, and all four links are removed, we 
have 



var 



£,[8-41] 



2Nt,Ni 



L iv R 



(A^L + A R ) 2 ( 2 P2-P3) 



With the tunnel barriers though, new possibilities also 
arise. For example, when a fringe tunnels into the lead, 
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one of the trajectories could be backscattered, as for the 
2-encounters previously. Similarly to (17), we obtain here 



var 



G [8-ll] 



8iV 2 iV 2 (l-r)p2 
(N L + N R ) 4 (2p 2 -p 3 ) 



x ( e -P3T E /r D _ e -2p 2 T E /r D j _ ^ 

Furthermore we can move additionally one fringe into the 
lead yielding 



var 



G [8-31] 



4(N h Nl + N£N R ) (l-F)F 2 
(N L + N R f (2p 2 -p 3 ) 



( 



>-P3T E It D _ p -2p2T B /T D 



(58) 



If both fringes are partially backscattered at the leads, 
we have the contribution 



var 



£.[8-21] 



8/V 2 N 2 

OJV L JV R 



r(i-r) i 



(JVl + JV r ) 4 (2p 2 -p 3 ) 



x I e 



-P3T e It d 



-2p 2 T E /T D 



). (59) 

Next we can allow the periodic orbit itself, rather than 
the fringes, to touch the lead and be reflected from the 
tunnel barriers. The encounter stretches may then both 
tunnel into the lead, or one may be reflected. We recall 
that the case in which both are reflected is already in- 
cluded in the survival probability of the encounter. One 
fringe is therefore removed from the systems, while the 
fringe that remains could then be inside the system, tun- 
nel into the lead or be partially backreflcctcd at the open- 
ing. If the periodic orbit touches both leads during the 
encounter, we return to the situation in the previous sub- 
section, so here we only consider the cases in which at 
least some of the second fringe remains. 

The removal of one fringe entirely implies the follow- 
ing changes in (53): With only one fringe time left, we 
only have one of the factors in the third line of (53). 
Also the encounter time is replaced by an integration 
variable (with a range from to t CDC ) which accounts 
for the amount of the encounter left after the periodic 
orbit is reflected from the tunnel barriers. Performing, 
with these changes, again the steps in [6] analogous to 
those following (53) above, we obtain the contributions 
for the various different cases. First, when the remaining 
fringe lies inside the system, we could have the encounter 
(which touches the lead) partially reflected 



varGt 8 -"- 11 ] 



(JV L + N K f P3 
(1 



P3TeIt D \ -p 2 T E /T D 



(60) 



where '-lf-11' refers to the fact that one fringe and one 
link were removed. Likewise, both encounter stretches 
could tunnel directly into the lead so that two links are 
removed on the same side of the encounter 



(N L N 2 + NlN R )(l-T)p 2 r 
(N L + N R ) 3 P3 
x(l- c - p3TE/TD )e- p2TE/TD . (61) 



The remaining fringe could also be paritally reflected 
from the tunnel barriers. For the encounter stretches 
have the same two possibilities as above. This leads to 

varG [8-if-2i] _ ^NlNl (i-r) 3 r 

(iVL + iVR,) 4 P3 

x(l-e~ P3rB/TD )e~ P2rB/TD , (62) 

and 



var 



G[ 8-if-3i ( s)] _ 2(^jv R + jv L jv R )(i-r) 2 r 2 

(JV L + N R f P3 
x(l- e - p3TE/TD )c- p2TE/TD . (63) 

Finally the remaining fringe can tunnel directly into the 
lead, giving the contributions 

varGt 8 - 1 ''- 31 ] - 2 (^R + ^LiV R )(l-r) 2 r 2 
(JVl + N R f P3 
x ( 1 - e - p3TE/TD ) e - p2TE ' TD , (64) 



and 



varG [8- lf -4i] = 2A^ R (i-r) 2 r3 

(N L + N R ) 2 P3 



x ( 1 _ e -p3TE/T D j e -p 2 TE/T D (g^ 



E. Multiple periodic orbit traversals 

Having treated all encounter configurations with pe- 
riodic orbits where the encounter stretches are shorter 
than the enclosed periodic orbit we now turn to the cor- 
responding contributions where the encounter stretches 
are longer than the enclosed periodic orbit. In this con- 
text only diagrams with stretches that both possess the 
same number of traversals around the enclosed periodic 
orbit yield a contribution, for a justification see [6]. In 
this context we consider that each encounter stretch has 
k full windings around the enclosed periodic orbit. We 
take here the two encounter times i cnc ,i, i C nc.2 and the 
corresponding primitive times t p nc i = t cnCt i - kr p . These 
primitive times are again shorter than the enclosed peri- 
odic orbit and we can now also consider the different cases 
depicted in Figs. 5-8, and how they can be arranged. We 
start with the case in which the two primitive encounter 
times, t p encil do not overlap. When k = 0, the corre- 
sponding diagram is depicted in Fig. 5. To explain the 
calculation we begin with Eq. (28) for k = and perform 
the integrals over the links U , and over t' where the latter 
integral leads to the factor (t p - t cncs ): 



varG 



[5] 
fc=0 



L R , / dsdM-| 



e* f 



(A^l + N R ) 

r- oo 

X / dT p (T p 

•J tones 



enc, 1 ^cnc,2 



^encs) 



XC 



encs J 

tones 

one;- 



(66) 
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To turn to k + 0, i nc,i is then replaced only inside 
the integrand of the r p -integral by t p nc v these times are 
shorter than t p by definition. The t cnc ,i before the t p - 
integral compensate the overcounting of equivalent po- 
sitions of the Poincare surface of sections by the Sj,Wj- 
integrals and thus do not need to be altered. The pos- 
sibilities for placing these stretches around the periodic 
orbit are the same as they were for i e nc,i for k = 0. This 
allows us to treat k > in essentially the same way we 



treated k = before. As t p 



= t„ 



kr p , we reex- 



press the primitive encounter times in terms of t enc ^ and 
T p . The limits of the r p -integration are also altered for 
k + 0: As the primitive encounters do not (yet) overlap, 
we have the condition t p enc 1 + t p cnc 2 > t p so that the lower 
limit is now t encs /(2k +1). Since we remove all k com- 
plete windings of the periodic orbit from both encounter 
times, the shorter encounter must be t C nc,min ^ kr p so 
that the upper limit is i C nc,min/&- Finally pj is replaced 
by P2k+j (where in this context T is denned as pi) since 
during the primitive encounter stretches we have a to- 
tal of 2k + 2 correlated stretches while elsewhere we have 
2k +1. All these replacements finally yield 



varG 



[5] 



NlN? 



J„\ 4 J-c 



(A L + JV R ) 
*[ 



dsdit 



Ti 



H 



e" 1 



n 2 1, 



enc,l^enc,2 



(67) 



dT p [(2k+l)T p -t cncs ] 

-[(P2fc + 2-P2fc + l)(tenc S -2fcTp)+p2fc+lT p ]/r£> 



xe 



where the additional index 'mt' indicates that this contri- 
bution results from multiple traversals k > 0. Performing 
now the remaining integrals in the same way as for exam- 
ple in Eqs. (34-38) we obtain this contribution. As the 
full expression is rather involved we just give here terms 
up to linear order in te- 



var 



G [5] 



2NINI 



_ y. (P2fc+2 -P2k+l)T E 

~(N L + N R fk fc 2 (2fc+l)r D 
+ 0(t%). (68) 



We also omitted here terms that are also nonzero for 
Pk = 1 as they all cancel with the contributions below. 
Also in this case one link can be removed leading to a 
configuration where one encounter stretch touches the 
opening. This implies that in (67) one factor td/T is 
removed, the duration of the encounter stretch which 
touches the tunnel barrier is replaced by an integration 
variable integrated from to t C nc.i 7 and additionally a 
factor (1 - Y) 2k+1 arises to account for the probability 
that all the 2k + 1 orbital parts surrounding the periodic 
orbit are backreflected when they hit the lead. This con- 
tribution then yields 



var 



r [5-ll] 
"mt 



4iVfM f 
(iV L + N R ) 4 ti 



(P 



2k+2 



■P2k+l)TE 



k 2 (2k+l)T D 



rl). 



(69) 



Shrinking two links connecting two different encounter 
stretches to the opening, we obtain (following the same 
steps as just described also for the other t e nc,i) no con- 
tribution which is linear in te- Shrinking two links con- 
necting the same encounter we need again to introduce 
in (67) two additional time integrals as described before 
(47). This also finally yields zero contribution linear in 
te, because the terms resulting from the two limits of the 
Tp-integral cancel. The same also holds for contributions 
following. 

Similar contributions are also obtained in the other 
cases. For a generalized 3-encounter introduced around 
Eq. (29) we obtain the integrals given in (Al) which pro- 
vide the contribution 



varG 



[6] 



2N t N K 



TE 



(N L + N R f ti Hk+1)T D 

X (2p 2 fe+2 ~P2k+l -P2fe+3) +0{t%) 



(70) 



when the encounters lie inside the system and 



var G 



[6-11] 



* N l N R 



TE 



(N L + N R f fct! k(k + 1)t d 

X (2p2k+2 ~ P2k+1 ~ P2k+Z) + O (r|) 



(71) 



when one link is removed. Additionally we can shrink 
both links connecting the 3-encounter to the opening on 
one side. This amounts to removing two links, and in- 
troducing an additional integral over the part of the en- 
counter that remains inside the system, in (Al), yielding 



varG 



[6-21] _ 



2(N L N 2 +N 2 N R ) 
(N L + N R f 

X (2p2k+2 ~P2k+l 



y IE ( 72) 

t,Kk+l)TD K ' 
P2k + 3)+0(T 2 E ). 



Again zero contribution is obtained when more links are 
removed. 

For a non-aligned 3-encounter, when the encounter lies 
inside the system we obtain from Eqs. (A2,A3) 



varG 



[7(a)] 



(A^l + A R ) 4 hi 



2fc+3 



■P2k+2)T E 



(k+ l)(2fc+ 1)t d 
(P2k+2 ~ P2k+l) T E 



k(2k+l)T D 
and when one link is removed: 



varG{I (a) - 11] 



8NIN 2 - 
(N L + N R ) 4 ti 



+ 0(t 2 e ) (73) 



(P2k+3 ~P2k+2)T E 
(k+l)(2k+l)T D 



(2p2fc+2 ~ P2k+1 ~ P2k+3) TE 

2/,(/,-+ \ )r n 

(p2k+2 ~P2k+l) TE 



k(2k+l)T D 



+ 0{t%). (74) 



Finally we treat the configuration in which the encounter 
stretches overlap at both ends. When k = this configu- 
ration is depicted in Fig. 7(b). With multiple traversals 
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of the periodic orbit, the contribution given in Eq. (A4) 
provides, when the stretches lie inside the system, 



varG 



[7(b)] _ 



(iv L + N R y 



E 

fe=i 



(P2fc+3 -P2k+2)l~E 

(/c+l) 2 (2fc+l)r D 



(75) 



and when one link is removed 



varG 



[7(b)-ll] 



4iV 2 iV 2 



E 



(N L + N R f fe tl(fc + l) 2 (2fc+l)r D 
+ 0(ri). (76) 



No further contributions are obtained when taking into 
account fringes as in Fig. 8. This can be checked from 
Eq. (53) by adjusting the limits of the r p -integral ap- 
propriately. Then, the terms from the two limits of the 
Tp-integral cancel to linear order in te ■ 



F. Linear Ehrenfest-time dependence 

Having calculated all the contributions we can now 
sum them to obtain the overall contributions to the con- 
ductance variance, to leading order in the number of open 
channels in the leads. First we can check that all the 
terms after Sec. Ill A are zero for te = 0, so that the 
RMT result (25) is preserved. The case that is especially 
interesting for the comparison with the numerics is the 
contribution to the conductance variance when Al = JVr, 
for arbitrary T and of linear order in te- We therefore 
consider for Al = Ar -»■ oo, in the presence of time rever- 
sal symmetry, 



te 



var G iv L = iv R ~ var G +a — 

td 



i[l + (l-r) 2 ] + a 



TE 
TD 



(77) 



We obtain 



a (r) = -[(r-i)(7r 3 -6r 2 + 4r-2) + in(2r-r 2 ) 



r 2 (2-r) 
: (i-r) 



-rLi 2 ((i-r) 2 ) 



(78) 



by summing the contributions linear in te from Eqs. (16- 
24, 26, 39-40, 43-45, 48-52, 54-65, 68-76). Here Li 2 (x) 
denotes the polylogarithmic function, which arises when 
performing the fc-summations above. The function a(F) 
is shown as the solid (red) line in Fig. 2. 

Having derived the complete set of contributions, we 
now discuss the main terms that lead to the two ex- 
trema in Fig. 2. First we note that the contributions 
from multiply-traversed periodic orbits are quite small. 
When T = 1 they provide no contribution linear in t e , as 
can be checked from Eqs. (68-76). For most of the con- 
tributions for r + 1 the index fc of the pk 's in the sum is 



much larger than one. In this case the differences of the 
Pfe's in Eqs. (68-76) tend to one and thus the contribu- 
tions from these equations tend to zero. Then the result 
is almost equal to the vanishing value for T = 1. 

In Fig. 2, the maximum lies in a region of small V. 
In this regime, most of the diagrams approximately can- 
cel each other as can be seen by expanding their con- 
tributions in a Taylor series for small T: The ones from 
diagrams containing encounters inside the system [e.g. 
the contribution (16)] are canceled by the contributions 
obtained when one and two links are shrunk [e.g. the 
contributions (17-21)]. However, this cancellation does 
not hold for the 3-encounter with the periodic orbit (see 
Fig. 6). In this case we have two possibilities: First that 
the encounter lies inside the system or touches the lead 
at one end [see Eqs. (22-24)]. In total this leads to an 
Ehrenfest-time dependent factor oi e~ P3TE ^ TD . Addition- 
ally, we also need to take into account the contributions 
from encounters touching the leads on both sides leading 
to an increasing contribution proportional to 1— e~ P3TE l TD 
with a larger prefactor than the first one [see Eqs. (50- 
52)]. Together these cause the first peak in Fig. 2. 

In the case of the dip the factor (1 -T), i.e. the proba- 
bility for the particle to be backreflected at the opening, 
is quite small. This implies that we get the main contri- 
butions from diagrams with encounters inside the system 
(i.e. the ones also obtained for T = 1) depicted in Figs. 
4-8. These together yield a negative contribution. 



IV. NUMERICAL SIMULATION 

We numerically confirm our main analytical prediction, 
Eqs. (77) and (78). The model we use is the open kicked 
rotator with time-dependent Hamiltonian [37] 



H = 



(P + Po) 



+ Kcos(x + x ) E ^( t ~ nT f)' (79) 



with Tf the free flight time. Depending on the kicking 
strength K the dynamics changes from integrable for K = 
to fully chaotic for K > 7. In the latter regime, local 
exponential instability is characterized by the Lyapunov 
exponent 



A: 



Tf^d) 



(80) 



The quantities p and x are introduced to break the 
Hamiltonian's two symmetries [37], and investigate dif- 
ferent symmetry classes with (p or x = 0) and without 
(po + + .To) time reversal symmetry The Ehrenfest 
time in this system is determined, within a constant of 
little relevance, by te = A -1 ln(M/c) with M the size of 
the Hilbert space, determined by the quantization of the 
Hamiltonian via discretization of the coordinates as, e.g. 
Pi = 2nl/M, 1 = 1,... M. Here c is a system-dependent 
constant of order one that is of classical origin, and as 
such does not depend on M. A quantum representation 
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varG 0.2 




0.15 



FIG. 9: Conductance variance var G(E) as a function of 
t~e/Ftd for various values of T. The straight lines are lin- 
ear fits constrained to go to (1 + (1 - T) 2 )/8 at te = 0. Their 
slope give the parameter a defined in Eqs. (77) and (78). 



Further dividing this slope by Y gives us a. There is 
still an uncertainty in a due to an uncertainty in the pre- 
cise value of the Lyapunov exponent - the latter has been 
found numerically to deviate in open systems from its ex- 
act value (80), which has been attributed to finite-time 
effects - and because of the arbitrary constant c of order 
one in the definition of te- We remove this uncertainty by 
forcing numerical and analytical data to agree for V = 0.2. 
Once this is done, there is no free parameter left. We 
compare the so obtained numerical values for a with the 
analytical prediction of Eq. (78) in Fig. 2. The excel- 
lent agreement between numerical data and the analyti- 
cal curve fully confirms our theory. We also checked, but 
do not show, that the linear r^-dependence of varG(E') 
is halved when time-reversal symmetry is broken. 



CONCLUSIONS 



of the Hamiltonian (79) is then provided by the unitary 
M x M Floquet matrix U , giving the time evolution for 
one iteration of the map defined by H in a time interval 
[to, to + Tf[. For our specific choice of the kicked rotator, 
the Floquet operator has matrix elements 

Uiv = Me' {m/Mm+lo)2+{l ' +lo)2] (81) 

E2irim(l-l')/M -{iMK/A-k) cos(2tt (m+m ) / M) 
e e 

m 

with Iq = poM /2n and mo = xqM /2tt. Transport can fi- 
nally be investigated once absorbing phase-space strips 
are introduced to model contacts to leads. This is 
achieved by means of projection operators P. In our 
case of tunnel-coupled leads, the latter are N x M diago- 
nal matrices with entries = Sij^/T, assuming that the 
system is coupled to all N = N-^ + Ar channels with the 
same transparency < T < 1. The scattering matrix is 
finally defined as [38] 

S(e) = (1 - P T pfl 2 - P[exp(-ie) -U(l- P T P)y 1 UP T . 

The kicked rotator model is particularly performant to 
investigate the semiclassical limit with te i5 To, as it 
allows for a rather large variation of the system size and 
hence the Ehrenfest time. Previous such investigations of 
the open kicked rotator are presented e.g. in Rcfs. [7, 9, 
35, 36], and we refer the reader to these works for further 
details of the model. The variance of the conductance is 
calculated by varying the quasienergies s at given lead 
position. Averages of var G(E) are further performed 
with different lead positions. 

We first show in Fig. 9 the behavior of var G(E) as a 
function of te/tjj. We clearly see T-dependent behaviors, 
as varG(iS) increases for small T and decreases for larger 
T until it becomes independent of te for T = 1. We ex- 
tract the linear slope of the varG(E') vs. te/td curve via 
a linear fit, setting te = where the curves intersect the 
RMT universal value var(G) RMT = (1 + (1 - r) 2 )/8 [28]. 



In this article we determined the dependence of the 
conductance variance of a chaotic cavity with tunnel bar- 
riers on the tunneling probability T. In particular we find 
an Ehrenfest-time dependence for the general case T < 1 
of non-perfect coupling. We focused on the contribu- 
tion which is linear in the Ehrenfest time and at leading 
order in the total number of open channels N. We pre- 
dict a nonmonotonous sinusoidal behavior with one max- 
imum at r Rf 0.2 and a (roughly 50%) deeper minimum at 
F «a 0.8. This analytical result was derived semiclassically 
by systematically considering a rather large number of 
possible configurations, but the general behavior derives 
from a much smaller set. We discussed the main contri- 
butions which lead to each of the two peaks to provide 
a better intuitive understanding of the overall structure. 
Finally, we compared these analytical predictions with 
numerical simulations performed for the kicked rotor (in 
the chaotic regime). There was good agreement between 
the two, with the analytical curve within the error bars 
of all the numerically determined datapoints. 

Although, even at leading order, we treated a large 
number of possible semiclassical diagrams, they could be 
constructed, as we showed in this article, in a controlled 
way from the diagrams that exist when the tunnel barri- 
ers are absent. This calculation could therefore possibly 
be extended to higher order contributions in 1/A and 
te, but a more natural next step would be to extend the 
Ehrenfest-time dependence of the (leading order in A) 
results for the full counting statistics of transport mo- 
ments [17] to include tunnel barriers in the leads and 
to account for the extra diagrammatic possibilities that 
then arise. 

On the numerical side, we emphasize that the agree- 
ment observed in this article is the first observed for the 
conductance properties of systems with tunnel barriers. 
Neither the results for the conductance nor for shot noise 
obtained in [34] could yet be confirmed by numerical sim- 
ulations. Further numerical investigation in this area 
would therefore be highly desirable. Finally it would be 
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interesting to check the observed phenomena experimen- 
tally. Although Ehrcnfcst-timc effects in antidot super- 
lattices were observed more than ten years ago [4], none 
of the other Ehrcnfcst-time dependencies predicted for 
chaotic systems could be checked so far experimentally. 
This has been mainly due to the small range through 
which the system size can be varied in quantum dots. 
However, with the additional parameter of the tunnel- 
ing probability and in view of the expected double peak 
structure, the prediction in this article could be amenable 
to experimental verification. 
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and second for <f nCjmax < t p < ^ nc l + ^ nc 2 
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2t d 
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x [e 



xe 
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Appendix A: Periodic orbit encounters 



In the case of the two encounters overlapping at both 
ends the contribution is obtained to be 



Here we give the necessary integrals for obtaining the 
contributions from diagrams containing enclosed periodic 
orbits that are surrounded fc-times by each encounter 
stretch. In these expressions the link and ^'-integrals in 
Eq. (28) have already been performed, but not yet the 
t p - and the s, it-integrals. These expressions are espe- 
cially useful because they not only allow one to obtain 
the contributions from encounters inside the system but 
(via integrations with respect to the corresponding en- 
counter times) also the contributions when encounters 
touch the openings. 

The corresponding expression obtained for two inde- 
pendent 2-encounters was already given in the main text 
in Eq. (67). In the case of a generalized 3-encounter we 
obtain 



var 
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L R . / dsdu I dr„-^ 



e^ 1 



(t. 



xe 



xe 



^cnc,l^enc,2 

(P2fc + 3-P2)s + 2)(t 







-[(p2k+2-p2k+l)(t. 



cnc,max L cnc,min ) 

kr p )/r D (A-l) 

kT p )+p 2 k + lTp]/T D 
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cnc,max 



For a non-aligned 3-encounter, similar to (34,37), it is 
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encs 
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Appendix B: Different tunneling probabilities 

Here we generalize Eq. (25) to the case of different Tj 
for the different lead modes. Factors like FN are then 
replaced by a sum over the Tj with respect to the N 
open channels. An analogous replacement is made for 
the pj for j > 1. In order to keep the notation compact 
we introduce 

Ni N L +N R , 

G (4) E^r„ Gi = Y, Vr (Bl) 
j=i j=i 

For example, the contribution from the diagram in Fig. 
4(a), previously given in (16), becomes 



var 



G [4] 



r<2 r 1 ^ r<2 



(B2) 
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where the Ehrcnfcst-time dependence is the same as in 
(16) although we set te = here. The first two terms 
in the numerator result from the channel summations. 
The third term replaces the factor p\N 2 in (16) and the 
denominator replaces the previous (IW) 4 . 

To keep the contributions from configurations where 
one link connecting the encounter to the opening is re- 
moved in a compact form, we define 



3=1 



fc-1 



(B3) 



Considering again the corresponding contribution from 
Fig. 4(a), previously given in (17), we have 



var( 



[4.U] 2 [ff(i), 2 G(2) + G (1 )g(2), 2 ] G 2 G(i)G( 2 ) 

= Gl • 

(B4) 

The first term derives from when the link connecting the 
encounter to lead 1 is shrunk while the second term corre- 
sponds to when the link connecting an encounter to lead 
2 is removed. Compared to (B2), the factor H^^ k takes 
into account that the particle is entering the system in 
a certain channel and returns to the same channel after 
traversing a link. 

Removing two links from different encounters likewise 
leads to 



var 



G [4(o)-21] = 4G (l) g (l),2g(2), 2 G(2) 

G\ 



For the diagram in Fig. 4(a) we can also remove two links 
from the same encounter, and for this we define 



Ni 



/(o^Er^i-r,-) 



fc-i 



(B6) 



3=1 



so that we obtain 



wgiw)] = [h^Gi 2) + G y {2) ,] G ^ (b7) 

G 1 

r ,[4(a)-3l] _ 2 [ G (l) H (l),2l(2),l + /(l),lg(2), 2 G(2)] 

var or - ^ , 

(B8) 

G [4(a)-41] = / (1)4 / (2),1 



var 



G\ 



(B9) 



For the diagram in Fig. 4(b) we cannot shrink two links 
attached to the same encounter (or more than two links) 
and because of the way the encounters are arranged when 
we shrink two links we obtain 



var 



G [4(6)-21] = 2G (l) g (l),2g(2),2G(2) 

G\ 



fr2 r~t2 , rt2 tt2 
J3 (1),2 Lt (2) + Lt (1)'"(2),2 



G\ 



Using the definitions above, this is equal to (B5) plus the 
results from (B7 B9) so that time reversal symmetry still 
leads simply to a factor of 2. 

In an analogous manner, we can show that each dia- 
gram in Fig. 5 gives twice the contribution of the diagram 
in Fig. 4(b) while for the diagrams in Fig. 6 we obtain 



var 
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2G (1) G (2) G 3 



G\ 
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var 



var 



r [6-U] _ 4 [ G (l) g (l),3 G (2) + G (l) g (2),3G( 2 )] 

= G~\ ' 

(B12) 

G [6-21] = 2 [ J (l) ,2 G (2) +g (l) / (2),2] 



G\ 



(B13) 



Summing all these contributions, we obtain the RMT 
result [28] for the leading order in TV contribution to the 
conductance variance 

var G RMT = 1 [2g\g' 2 - 2 g \g' 2 g\ - 4gfg' 2 g[ 2 
-45iff25i -2ffig 2 5i + 2 92ffi +2g 1 g 1 

, A 3 '3 „ 2 ,4 „ 4 ,2 - 2 , ,2 

+ 4 5i5i + 2 5i3i +35i52 +»9i92g 2 9i 
+ 3glg[ 4 - 2g\g' 3 - 2g i 1 g' 3 g' 1 - ^gig^ 

- 2 535i 5 ] , (B14) 

for systems without time reversal symmetry, and twice 
this result for those with. Here the notation introduced 
in [28] 



was used. 



9k = Y, F "> 9k = E r n> 

n=l n=l 



Appendix C: Shot noise 



(B15) 



Here we calculate the shot noise power. It can be writ- 
ten as [39] 

P = G-h= (Tr[ttt])-(Tr[(ttt) 2 ]). (CI) 

The average conductance G(E) has previously been cal- 
culated [34, 40], and the first few terms for systems with 
time reversal symmetry can be written as 

G{E) __ G^^^G^G, (C2) 

(2 \ G(l)g(2),2 +g(l), 2 G(2) 

1/3 ) G\ 

using the notation of Appendix B. The second two terms 
derive from the diagrams in Fig. 3 so that without time 
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FIG. 10: Two trajectories with a single encounter and two 
partner trajectories. 



reversal symmetry only the first term in (C2) remains. 
This is included as j3 = 2 for systems without time rever- 
sal symmetry and j3 = 1 for those with time reversal sym- 
metry. The function h is given semiclassically in terms 
of 4 trajectories 



iE E E A 1 A* 1 ,A i A\, 

1 H a,b ■y(a^b) {(c^rf) 
c,d 7 '( c ^fc)£'(a^d) 



xe 



■k( s -,-S- r ,+s i -s i 



■">■ 



(C3) 



where the main difference from the conductance variance 
(7) is that the trajectories connect different channels and 
the sum is unrestricted. We start with the diagram from 
Fig. 10, where since the partner trajectories cross over in 
the encounter they automatically travel from and to the 
correct channels and all four channels are unrestricted. 
This means that this structure is now lower order (in in- 
verse channel number) than it was for the conductance 
variance (it was there contained by taking into account 
the diagram in Fig. 4(a) and removing two links con- 
nected to the same encounter). For the diagram in Fig. 
10, we therefore have 



G2 ri 2/~t 
(1) ( - T (2) ^2 



2 [ G (i) g (i),2 G ( 2) + G 2 (1) H (2)a G (2) j 
G\ ■ 



(C4) 



This contribution, along with its Ehrenfest-time depen- 
dence, was previously calculated in [34]. 

Moving to the next order term, which has recently been 
calculated using RMT [41], we again simply need to look 
at the corresponding structures and consider the possible 
ways of shrinking links. Without time reversal symmetry 
there are no possible structures at this order, so we there- 
fore consider systems with time reversal symmetry. The 



structures which contribute are depicted in Fig. 2 of [15]. 
There is also, in Fig. 2(d) there, a quadruplet involving 
two independent pairs, one of which is simply a diagonal 
pair, while the other involves a single 2-encounter, as in 
Fig. 3 of this paper. Because the start (or end) channels 
must coincide a = c (or b = d) there is a subtlety when we 
consider shrinking one link. If we shrink the link on the 
left (or right) side of the 2-encounter so that it moves into 
the lead where the diagonal pair emanate (or terminate) 
we actually have the same case as when the 3-encounter 
from Fig. 2(c) (of [15]) moves into the lead. As such, we 
include this case there, leaving a contribution of 



kPCd)] 



! [ G (1) 7 (2),1 + i 



G (2)] G 2 



G\ 



(C5) 



2 [G , ( 1 )ff( 1 ) i2 /(2),l + 1(1), 1#(2), 2 G (2)] 



G l 



The diagram in Fig. 2(a) of [15] is similar to the diagram 
in Fig. 4(b) here and we obtain 



4G (1) G (2) G 2 + 8G(i)-ff(i) i2 ff(2) : 2 G (2) 



G l 



G\ 



8 [G(1)#(1),2 G (2) + G (i)#(2),2 G (2)] G 2 



'■[ G h H f 



Gl 



2 , tt2 (~i2 
2), 2 + X1 (l),2^ r (2) 



Gf 



(C6) 



For the diagram in Fig. 2(b) of [15] we obtain twice the 
result of Fig. 5(a) here [or four times the result of Fig. 
4(b)] for the conductance variance. For the diagram in 
Fig. 2(c) of [15] we likewise obtain twice the result of Fig. 
6(a) here for the conductance variance. With all the next 
to leading order shot noise contributions, we can combine 
them and indeed find the same result as in [41]. 

We can simplify the result by setting all of the tunnel- 
ing probabilities equal to T: 



P(T) 



r(i - r)NiN 2 + r(3r - 2)NfN% 



(N h + N R ) 



V/3 



T{AT-2,)NiN 2 {Ni- N 2 ) 
(JVl + iV R ) 4 



(C7) 



Setting r = 1, we also recreate the first two terms, 
N 2 N 2 

p(r = i) = iVliV2 , (cs) 

(^Vl + ^r) 3 

12 A( NxNi{N x -Ni) 2 \ 
+ ){ (N L + N R ) 4 j + 

of the result in [15]. 
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